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Abstract 

We carry out a theoretical and Monte Carlo study on the J/ip 
decays into ojirir and KKtt through intermediate subthreshold p meson 
by using SU(3)-symmetric Lagrangian approach. It is found that the 
subthreshold p contribution is not negligible and may have significant 
influence on partial wave analysis of resonances in these channels, 
especially near the urn and KK thresholds. 



1 Introduction 

In recent years, much effort has been devoted to the study of meson and 
baryon spectra. Of particular interests, the J /if) decays now provide an 
excellent source of information for studying light hadron spectroscopy and 
searching for glueballs, hybrids, and exotic states [Tj. Several interesting 
near-threshold structures were observed and studied, such as the broad a 
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near nir threshold j2], the broad k near Kir threshold the narrow /o(980) 
peak near KK threshold jl], and the narrow structure near pp threshold 
[Sj. While various mechanisms were proposed to explain these structures 
the conventional t-channel meson exchange final state interaction mechanism 
[7J IB] can give consistent explanation for all these structures. The t-channel 
p meson exchange was found to play a very important role for all these 
structures. 

In this paper, we want to address two other puzzling near-threshold phe- 
nomena in J/ijj decays. The first one is the "bi puzzle" in the J/ip — * unit. 
Both DM2 Collaboration [Hj and BESII Collaboration j2] obtained a much 
broader width for the near-u;7r-threshold resonance &i(1235): while PDG 
gives the width (142 ±9) MeV QU|, the DM2 and BESII gave it as (210 ±19) 
MeV and (195 ± 20) MeV, respectively. The second puzzle is that there 
is a clear near-^fT-threshold enhancement in both J/ip — > K + K~ix° and 
J/ijj — > KsK ± tt t Dalitz plots from DM2 data ^T] although the structure 
was not addressed in the paper possibly due to uncertainty of background 
contribution. From conservation laws for strong interaction, the KK here 
should have isospin 1 and spin-parity 1~. There is no known resonance of 
these quantum numbers very close to KK threshold. 

It is rather tempting to claim some new near-threshold resonances here. 
But before claiming any new physics from the seemingly puzzling phenomena, 
one should investigate all possible conventional mechanisms to see if the 
phenomena can be interpreted within the existing theoretical framework. 
Motivated by this idea, here we investigate the subthreshold p contribution 
to these channels through diagrams shown in Fig^and FigJ2]for J/ip — > amir 
and J/i/j — ► KKtc, respectively. There are good reasons for considering this 
mechanism. The J/ip — > pn decay has the largest branching ratio among 
the known two-body decay channels of the J/ip [10J. The pKK and pum 
couplings are well determined to be large. 

This paper is organized as follows. In section II, we present the method 
and formulation for calculation and Monte Carlo simulation. The numerical 
results and discussions are given in section III. 

2 Method and formulation 

The Feynman diagrams for relevant processes with the subthreshold p con- 
tributions are depicted in FigJT] ( J/xjj — > ojh + -k~ ) and FigO (J/i/j — > KKn ). 
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Figure 1: Diagrams for J /if) — * ujtt+tt decay with subthreshold p exchange. 




Figure 2: Diagrams for J /if) — ► KKtc decay with subthreshold p exchange. 



For the vector- vector-pseudoscalar (VVP) and the pseudoscalar-pseudoscalar- 
vector (PPV) couplings, we use the SU(3)-symmetric Lagrangians as in [51 IT2] 

C-WP = —^{d^daVpP), (1) 

Upv = - l -iG'{[P,d»P)V»), (2) 

where (. . .) means SU(3) trace, G and G' are the coupling constants, and 
P is the 3x3 matrix representation of the pseudoscalar meson octet, here 

P = \ a P a , a = 1, ,8 and A a are the 3x3 generators of SU(3). A similar 

definition of V v is used for the vector meson octet. 

In the Gell-Mann representation, the relevant effective Lagrangians are 

/V = (3) 

C*„ = g^e^p^(u) P y(p), (4) 

Cpnir = g P nA(P«+ ~ + ~ P£o)p£ + (P^o ~ P^+)Pp] (5) 

£ P KK = g P KR[{PK+ ~ PK-)PI + (PRO ~ + 

g pK K{p K o - P K -)pl + v 7 ^ g pK K{p K + - Pro)p p ( 6 ) 

where <7 p #/f = ^p?™ = C' due to flavor SU(3) symmetry. Using these La- 
grangians, we are able to construct following amplitudes T\, T 2a , T 2 6, T 2c and 
T 2 d corresponding the diagrams in FigJTJ FigJ2k, FigJ2b, Fig|2fc and FigOl, 
respectively: 

Tl = -giipirgupirUtMvapPoe*' ' (if)) (pi + p 2 )° ''^^ '{j>\ + P2) 7 Pl\e a {v) 

1 1 

+ 



(pi + P2) 2 -m 2 p + impT p (pi + p 3 ) 2 - m 2 + 2m p r p 

e All ,Q/3Poe*' / (^)(pi + p 3 ) a e 7/3A<T (pi + Pa^lAe^)}^^^, (7) 

T 2a = r a = -^r te = -^r M 

2g<ppTrg p KK n *v( i\ a (3 



(Pi +P2) 2 -m 2 + im p T^ 



»vcLp^Z* V (^)V\ThP^(mFpKK- (8) 



Here T p represents the width of p; F^ p7r , F^p^ and -F pii -^- are the form factors 
for the ippTt, ujpn and pKK vertices, respectively. 

Usually, hadronic form factors should be applied to the meson-meson- 
meson vertices because of the inner quark-gluon structure of hadrons. It 
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is well known that form factors play an important role in many physics 
processes, such as tttt scattering J3J, NN interactions scattering 
|14j . meson photo-production [H^ etc.. Due to the difficulties in dealing 
with nonperturbative QCD hadron structure, the form factors are commonly 
adopted phenomenologically. 

The most commonly used form factors for meson-meson-meson vertices in 
J/ip decays are Blatt-Weisskopf barrier factors Bi(Q abc ) [T^J HZj for a decay 
process a —* b + c with orbital angular momentum I between b and c mesons: 

B 1 (Q abc ,R) = \ L 2 Q l n2 for/ = l, (9) 

V "abc + ^0 

where 

Qlbc = {Sa+ 2~ Sc)2 - s b (10) 

is the magnitude of pb or p c in the rest system of a; Si = Ef — with 
Ei and p; the energy and the tree- momentum component of pi, respectively. 
Here Q is a hadron "scale" parameter, Q = 0.197321/.R GeV/c with R 
reflecting the radius of the centrifugal barrier in fm. We take R = 0.5 fm for 
ippK vertex and R = 0.5 or 0.8 fm for upn and pirn vertices. For the sake 
of convenience, we use the following shorthand notation for the F^ pv F pram in 
Eqs. (j7l8|) for various channels: 

{B^Rf^ = 0.5)Bx(R fm ^ = 0.5) , J/ip -A tttttt 
B 1 (R^ fm = Q.h)B l {R W(m = Q.h) , J/^Awtttt (11) 
Bi(i2^,r = 0.5)fli(i^ = 0.5) , J/^^KKti 



BtiR^ = Q.h)B x {R^ = 0.8) , Jji) A TTTTTT 
B x {R il(m = O^B^R^ = 0.8) , J / ip —* ujtxtx (12) 
B 1 (R^ jp7r = 0.5)B 1 (R pKR = 0.8) , J/^j^KKtx 



The monopole form factor is also a frequently used s-channel form factor 

[El HUGH]: 

A 2 _|_ m 2 

FlA - q) = lsT7' (13) 

where m and q are the mass and the four-momentum of the intermediate 
particle, respectively, and A is the so-called cut-off momentum that can be 
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determined by fitting the experimental data. We take a commonly used value 
A=1.5 GeV for pirn [H] and upn ^H] vertices. For pKK vertex, A=1.5-4.5 
GeV was used in literature [7J |S] . It is possible that mesons involving heavier 
quarks have smaller size and corresponding need larger cut-off A parameter 
for the relevant vertices. We take A=l. 5-4.5 GeV for J /ippir vertex and pKK 
vertex. Similarly, we define 

[ F(A i>fm = 4.5)F(A p _ = 1.5) , J/V> A tttttt 

F 3 = I F(A^ = 4.5)F(A W/OT = 1.5) , J/ip-^ ujnn (14) 

[ F(A^ p7T = A.5)F(A pKR = 1.5 - 4.5) , J/V A KKn 



F(A^ = 1.5)F(A p7r7r = 1.5) , J/V> ^ TTTTTT 

F(A^ = 1.5)F(A W/W = 1.5) , J/iP^lutttt (15) 

F(A^ = 1.5)F(A p ^ = 1.5-4.5) , J/ip KKn 

The differential decay widths can be evaluated as 

dT(J/^ - CUTT+TT-) = ^-|Tl| 2 d$3(P0;Pl,P2,P3), (16) 



2Afy 



(2tt 



r'l 



dT(J/iJ> ^ KKir) = ^(\T 2a \ 2 + \T 2b f + \T 2c \' + \T, 



2M^ 



2d\ 



with being the mass of J/ip and the three body phase space factor 

f \ x4/ \ rf3 Pi rf3 ^2 d 3 p 3 / 10 x 

^3(P0; Pl ,P 2 ,P3) = « (P„ -„ - ft - Ps) ^ f2Ei (2?r) 3 2£ , 2 (2jrf2B3 . US) 

^From above equations and following the covariant tensor amplitude method 
described in detail in Refs.jTJj, the entire calculations are straightforward al- 
though tedious. The numerical results are given in the following section. 



3 Numerical results and discussions 

In the model described so far, there are four relevant coupling constants, 
Sty/™-) Quipn and g pn7T (see Eqs. (J3HBJ)). The g p7T7T can be obtained by evaluating 
the process p —>■ tttt with various form factors. In an analogous way, we can 
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obtain g^ pn through the sequential decays ifj — > pn with p — > txtx and g WfrK 
through u) — > pi\ with p — > 7r7r by using the similar approach described in the 
section above. The experimental data of p — > 7r7r, ^ — > p7r, co> — > 7T7T7r decay 
widths are from the PDG [TO]- The results are shown in Tabled 

We use two forms of T p in Eqs. ()7l8|) . First we take it as a constant of 
0.149 GeV. Then we take an energy dependent width including tttt, uir and 
K K channels 



Tp(s) — r p _ >7r7r (s) + r p_ >a)7r (s) + T p _ tKK (s) = T p _ t ^ riT+um+KI q, (19) 

where energy dependent partial width Tp^n, Tp^Tr and r p ^kk are obtained 
by the method similar to Eqs. (jl6ll7|) . Usually only p — > mi width is con- 
sidered in the p propagator. But in reality, when the invariant mass of the 
off-peak p meson goes above ujtt and KK thresholds, the off-peak p can also 
decay into these new channels and the corresponding partial decay widths 
should be included in its total energy dependent width. This will make the 
tail of p propagator drops faster for energies above the new thresholds than 
the usual propagator with constant width. We find that the two forms of 
T p have little influence on coupling constants. So we use the same coupling 
constants in amplitudes with the same form factor (see Table Q). 

In order to demonstrate how large effect the s-channel subthreshold p 
exchange may have upon various channels, their contributions are calculated 
and compared with experimental decay widths of the corresponding final 
states. The results are listed in Table [T] The R^ nn and R P K ^ represent 
the ratio of theoretical contribution from the subthreshold p exchange to the 
experimental width taken from PDG ^U] for the channels u)n + 'K~ and KKir, 
respectively: 

dp = T th (J/j; ^tt+tt-) 

The range of variations for R P K f Cn in Tabled for form factor F3 and F4 comes 
from the variation of A between 1.5 and 4.5 GeV for the pKK vertex. 
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Table 1: The coupling constants and the percentages of the s-channel sub- 
threshold p exchange contribution as functions of the form factors (F. F.) 
and of the choice of p width. 



F. F. 




Q pirn 


9,J/9pn 










(GeF) 




(GeV- 1 ) 


(GeV- 1 ) 


(%) 


(%) 


F x 


0.149 


8.18 


8.61 


14.56 


21.0 


4.5 








xl0~ 3 




18.3 


4.3 


F 2 


0.149 


10.69 


8.70 


14.68 


11.1 


4.2 








xl0~ 3 




9.3 


4.0 


F 3 


0.149 


6.05 


2.31 


11.82 


7.7 


1.7-3.4 








xl0~ 3 




6.3 


1.6-3.2 


F 4 


0.149 


6.05 


2.35 


11.82 


3.9 


1.0-1.8 








xl0~ 3 




3.1 


0.9-1.7 



From Table ^ we see that no matter which form of T p and form factor 
are employed, the range of the R^ n7T is from 3.1% to 21.0% and that of R P K ^ 
is from 0.9% to 4.5%. It means that the contribution of the s-channel sub- 
threshold p exchange is not negligible for both J /if) — > ujtttt and J/ if) — > KK-k 
channels. The subthreshold p contribution may have significant influence on 
the analysis of resonances near urn and KK thresholds. 

In order to see the influence it may have on analysis of resonances near 
thresholds, we perform a Monte Carlo simulation to give predictions on var- 
ious invariant mass spectra and Dalitz plots for these two channels as shown 
in FigOJl FigJU and Fig0 in which we take form factor F 2 and r p =0.149 
GeV. The dotted lines in the invariant mass spectra denote the uniform 
phase space distributions without considering the dynamical interactions. In 
FiglBl a clear enhancement about 1.2 GeV near the ujtt threshold appears in 
both invariant mass spectrum and Dalitz plot. Comparing with experimen- 
tal results (Fig. 7 and Fig. 10 in Ref.jH]; Fig.l and Fig. 2 in Ref.j2J), one can 
expect that this enhancement, as the background of &i(1235), should reduce 
the measured width of &i(1235) from this reaction and may well explain the 
"&i puzzle". 

One interesting phenomena in the Dalitz plot of Figgis that there is also 
a clear enhancement band near threshold, which is shown more clearly 
in the FigEJfor the ttti invariant-mass distribution of J/if> — > um + 7T~ decay 
(solid line). This enhancement band comes mainly from the interference 
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1 1.5 2 2.5 3 2 4 6 8 10 

M(w7T + ) (GeV) M 2 (om + ) (GeV) 



Figure 3: The uoir invariant-mass distribution (solid line) and the Dalitz 
plot for J /if) — > dJ7T + 7r~ decay through p exchange with form factor F<i and 
r p =0.149 GeV, compared with phase space distribution (dotted line). 



terms of two Feynman diagrams in Fig^ because if we do not include the 
interference term of Fig. Q the visible bump located in 0.5 — 1.0 GeV will 
be much reduced in the tt + tt~ invariant-mass distribution as shown by the 
dot-dashed line. Therefore, the subthreshold p contribution may even have 
influence upon the analysis of the a meson from this channel. 

In FigJSJ a clear peak around 1.1 GeV near the threshold is also seen in 
the KK invariant-mass distribution. By looking at Table and Fig|3 it is 
natural to expect that the subthreshold p contribution is an important source 
for the near-threshold enhancement found in J /if) — > K + K~tt° decay by the 
DM2 Collaboration [TT]. Hence the result provides another evidence of the 
important role played by the subthreshold p contribution. 

Since the mass of p is more than 150 MeV below the urn and KK thresh- 
olds and the width of p is not very broad, some people naively assume that 
the subthreshold p contribution can be neglected in these channels of J/if> 
decays. However, this paper should greatly change this point of view for 
the channels with final state particles coupling strongly to the subthreshold 
p. A similar important subthreshold contribution was previously noticed for 
J /if) — > NNw channels from subthreshold nucleon pole |20j . 

In summary, the p exchange plays a very important role in many low- 
energy strong interaction processes, such as tttt scattering, 7rK scattering, 
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500 - 




MOjVI (GeV) 



Figure 4: The nit invariant-mass distribution for J/ip — > ujtt + tt~ decay 
through p exchange with form factor F 2 and r p =0.149 GeV (solid line), 
compared with the phase space distribution (dotted line) and the result ig- 
noring the interference effect between subthreshold p + and p~ contributions 
(dot-dashed line). 
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1 1.5 _ 2 2.5 3 2 _ 4 6 

M(KK) (GeV) M 2 (Ktt) (GeV) 



Figure 5: The KK invariant-mass distribution (solid line) and the Dalitz 
plot for J/ij} — > KKti decay through p exchange with form factor F 2 and 
r p =0.149 GeV, compared with phase space distribution (dotted line). 



ttN interaction etc.. In this paper, we extend the mechanism to interpret 
some long standing problems observed in J/tp decays. It is found that the 
subthreshold p contribution is not negligible for both J/ijj — > ujtxtx and J/ip — > 
KKn channels, and should be included in analyzing these channels. It may 
well explain the longstanding b\ puzzle near um threshold in J /if) — > ujtttt 
and the near KK threshold enhancement in J/ip — > KKn decay. 
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